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Supernova remnants have long been regarded as sources of the Galactic cosmic 
rays up to petaelectronvolts, but convincing evidence is still lacking. In this work 
we explore the common origin of the subtle features of the cosmic ray spectra, 
such as the knee of cosmic ray spectra and the excesses of electron/positron fluxes 
recently observed by ATIC, H.E.S.S., Fermi-LAT and PAMELA. Numerical cal- 
culation shows that those features of cosmic ray spectra can be well reproduced in 
a scenario with e + e~ pair production by interactions between high energy cosmic 
rays and background photons in an environment similar to the young supernova 
remnant. The success of such a coherent explanation serves in turn as an evidence 
that at least a portion of cosmic rays might be accelerated at young supernova 
remnants. 
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Subject headings: cosmic rays — supernova remnants — acceleration of particles 



Ever since the discovery made by lKulikov Sz Kristiansenl (119581 ). the knee of the cosmic 
ray (CR) spectra observed ~ 4 PeV has remained a puzzle for half a century. The origin of 
the knee has become a key problem of CR physics, as it is closely related to the acceleration 
of the Galactic CRs. There have been extensive studies on the issue (e.g., see the review of 



Horandell (120061 ) and references therein), however, no consensus has been reached yet due to 
the lack of high precision experimental information. 

Recently, the ATIC expe riment measured t he cosmic electron (electron and positron) 
spectrum up to several TeV (IChang et al.l 120081 ) with high precision and showed evidence 
for an obvious excess of flux above tens of GeV. In addition, the ground based atmospheric 
Cerenkov telescope H.E.S.S. and the space satellite Fermi, which are dedicated to high 
energy 7-ray detection, also reported the i r resu lts on the electron spectrum measurements 
( lAbdo et al.l l2009i ; lAharonian et al.l I2QQ8I . |2009| ). Despite some discrepancies between the 
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results of those experiments, all of them seem to show the excess of electrons compared with 
the background estimation. At the same time, the PAMELA experiment also pub lished 
their results on positron fraction for energies up to ~ 100 GeV (lAdriani et al.ll2009l ). The 
PAMELA results show clear po sitron excess, in agreement with the anom aly observed before 
by balloon-borne experiments (IBarwick et al .11 1 9 9 71 : lAguilar et al.ll2007l ). 



It is su ggested that the rise of the positron fraction may require a primary source of 
e + e~ pairs (jSerpicol 120091 ) . As early as the 1990s, Aharonian and Atoyan proposed a model 



involving the e + e~ pair production mechanism from the interactions between high energy 
7-rays and optical or ultraviolet radiation in the vicinity of discrete s ources, to explain the 
electr on / positron excesses from early time balloon-borne experiments (lAharonian fc Atoyan 
199ll ). We think that it might be more probable to explain the lepton excess by the e + e~ pair 
production via interactions between CRs and background photons. Such a process would be 
much more efficient than the 77 interaction because the cross sections are comparable, while 
in the former case multiple collisions are allowed. 

In addition to the production of e + e~ pairs, the interaction between nuclei and back- 
ground radiation will unavoidably change the spectrum of CRs on the source and may lead 
to the formation of the knee. Considering the radiation at the source is at the optical energy 
level (~eV), we expect the pair production to occur at ~ 1 PeV for protons and a few PeV for 
Helium, which just corresponds to the "knee" of CR spectra. In the rest frame of CR nuclei, 
the secondary electrons/positrons have energies around MeV, which turn to be about TeV 
in the laboratory frame. It is striking to note this is just the energy range where the excesses 
have been observed by ATIC, Fermi and H.E.S.S.. As for the energy budget, the energy 
density of the excess electrons observed by ATIC is about 3 x 10" 5 eV cm -3 between 0.1 and 
1 TeV, which is of the order of CR energy loss assuming a spectral break from —2.7 to —3.1 
at energy ~ lPevH. It can be seen that the interaction between CR nuclei and background 
photons naturally bridges the knee of CR spectra and the excesses of electrons/positrons. 

It is generally believed that the Galactic CRs are accelerated by the shock waves in su- 
pernova remnants (SNRs). We will focus on the acceleration of CRs by the neutron star inside 
the shell of young SNRs for our purpose, because that is where abundant hig h energy CRs 



(proje c tile) and background ph o tons (t arget) can be found . As described in iGaisser et al. 



(119871 ); iBerezinskii fc Ginzburgi (119871 ); IGaisser et al.1 (119891 ). the CRs up to ultra high en- 
ergies (e.g., 10 18 eV) can be accelerated in a very short period (~yr or shorter) after the 



1 Note that this estimate is based on the observational spectra. The energy density of CRs at the source 
will be much larger than that needed of the excess electrons due to the fact that the diffusion time of PeV 
CRs is shorter than the energy loss time of TeV electrons in propagation. 
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explosion of the supernova. The radiation field has an initial temperature of about several 
thousand K, which is equivalent to a photon energy of ~eV. Before the radiation field decays 
and the temperature cools down in hundreds of days, the e + e~ pair production process can 
occur many times for one CR particle. As both the magnetic field and radiation field are very 
strong in the case of young SNRs, the electron energy can not exceed a few TeV because of 
the synchrotron depletion and the inverse Compton scattering. It is worth no ting that there 



should also be energy gain for electrons through acceleration. As shown in IVannoni et al. 



( 120091 ) the electron acceleration in radiation dominated environments should exhibit a pile 
up spectrum a round the cutoff energy, which should be a few TeV in the case of young SNRs 
( IGaisserl Il990l ). The properties of young SNRs seem to fit well with what is required to 
understand the knee of CR spectra and the abnormal spectra of electrons/positrons, but a 
quantitative study is necessary. 

For simplicity we decouple the acceleration and interaction processes in the calculation, 
i.e., the nuclei are first accelerated to very high energies and then they interact in the 
photon field. Finally, we inject the nuclei and electrons/positrons in the Galaxy to calculate 
the propagation effect and derive the observational spectra on the Earth. Monte-Carlo (MC) 
method is used to calculate the energy losses of the nuclei and the energies of the generated 
e + e~ pairs. 

There are three processes of interactions between CRs and photons, i.e., pair production, 
photodisintegration and photo- pion production. T he c ross sections 



'uget et al 



photodisintegration are given in iBlumenthall (119701 ) and 
pion production cross section for proton is adopted from 
an A 0,91 dependence for other nuclei with atomic number At lStanev et al.lll985l ). 
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(2008 ), an d we employ 



We choose a power-law spe ctrum with in dex about 2 for primary CRs, according to the 
Fermi acceleration mechanism (IGaisserl Il990l ) . The actual spectral index of the three main 
components, proton, Helium and Iron, are slightly tuned in the calculation to match the 
data. For other nuclei which do not play a significant role we adopt the parame terization 
'jz = a + bZ — 0.6 with the fitting parameters a = 2.69 and b = —2.07 x 10 -3 (IHorandel 



20031 ). Here the 0.6 subtraction from parameter a is done to explicitly take into account the 
propagation correction from the observations to the source spectra (see below). The relative 
abundances o f various kinds o f nuclei after propagation are normalized to the fluxes at 1 
TeV given in iHorandell (120031 ). The propagation effects of nuclei are approached using a 
leaky-box model. The escape time for CR nuclei is taken as t csc (R) ~2x 10 s ° 6 y r ; 
where R = p/Ze is the rigidity of nuclei. This relation is fitted from the low energy B/C 
data (IPutze et al .1120091 ) and we extrapolate it directly to the PeV energy range. 



As for the radiation field around the source, the temperature evolution with time is a 
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crucial factor to define the dominant type of interactions between CRs and background pho- 
tons. In a simple case, the temperature remains unchanged when CRs are being accelerated 
to ultra high energies. The three types of interactions may then happen simultaneously. 
However this may not be true providing that the threshold energies of photodisintegration 
and pion production are much higher than that of pair production. The pair production 
should start earlier and slows down the acceleration of CRs, which makes the other two 
processes with higher thresholds more difficult to occur. The observational data actually 
support this argument. As shown by the dashed lines in the upper two panels of Fig{TJ if we 
allow the photodisintegration to be active when fitting the Helium spectrum, we find that 
protons are over produced and the calculated proton spectrum does not agree with the data. 
The suppression of photodisintegration c an be understood by considering that the radiation 



field decays and/o r cools down rapidly (IGaisser et al.l 119871 ; iBerezinskii fc Ginzburgl 11987 



Gaisser et al.lll989l ) compared to the time scale of the CR acceleration. This leaves only the 
pair production as the main interaction between CRs and radiation photons. Therefore in 
the following calculation we forbid the photodisintegration and pion production and only 
consider pair production. 

Based on the previous considerations, the temperature and effective photon density 
are fitted to the observed spectra of individual elements. As shown by the solid lines in 
Fig{TJ the calculated spectra of proton, Helium, Iron and the all particle CRs agree well 
with the measurements. In Table [I], the parameters used in the simulation are summarized: 
the temperature, the effective column density of photon and interaction time assumed in 
a blackbody radiation field, and the initial power law index of each element. In this work 
we assume the second knee (~ 300PeV) of CR spectra is due to the interaction between 
Iron and background photons but this assumption is not a mandatory one (see below the 
discussion on systematic check). The temperatures and photon densities of other nuclei are 
derived using a linear interpolation with respect to the atomic number A between Helium 
and Iron. The heavy nuclei (low charge over mass ratio) correspond to low temperatures 
because they need longer time to be accelerated to reach the pair production threshold when 
the temperature of the source has decreased. It can also be seen from Table [1] that if the 
photon field is as intense as a blackbody radiation, an interaction time from a few weeks to 
a few months is enough to produce the observed spectrum distortion. 

The generated spectra of electrons/positrons at the source location for individual nu- 
clei and all CRs are shown in FigfJl The propagation effect of electrons/positrons in the 
Galaxy has to be considered to determine the spectrum on Earth. Unlike CRs, the dominant 
effect for electron propagation is the energy losses due to inverse Compton scattering and 
synchrotron radiation at ene rgies higher than ~ lOGeV. In this work the GALPROP code 



fjStrong fc Moskalenkdll998l ) is used to calculate the propagation of electrons and positrons. 
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As shown by the lower thick lines in Figj2], the propagated electron spectrum becomes 
softer due to the cooling effects from en ergy losses. It is ho wever still harder than the ex- 
perimental data. But as pointed out in iGaisser et al.l (119891 ). the high energy tail (>TeV) 
of electrons could be depleted in the environment of young SNR by the strong magnetic 
and radiation fields. An accurate calculation needs to involve a detailed modeling of radia- 
tion, interaction and acceleration. As this work focuses on the interaction rather than the 
acceleration and radiation of electrons, we will simply adopt an exponential cutoff on the 
injection spectr a of electrons /positron s with a value experimentally suggested by H.E.S.S. 
measurements (lAharonian et al.ll2008l ). The low curves in Figj2]show the effects of different 
values of cutoff energy on the electron spectrum. 

The calculated (e + +e~) spectrum and the positron fraction for a cutoff energy E c ps 2 
TeV are shown by the thick lines in Fig|3] together with the observational data. The calcu- 
lated flux of (e + +e~) is scaled to match the data by a free normalization factor (quantita- 
tively, it is ~ 3% for this model configuration). This parameter represents the dissipation 
and escape probability of electrons/positrons in the source region. In Figj3l the term "bkg" 
represents the standard background contribution of primary electrons and secondary elec- 
trons/positrons from CR interactions with the interstellar medium. For t he background 
calculation, we use the diffusion + convection (DC) model (jYin et al.ll2009i ) with parame- 
ters tuned to be consistent with the low energy CR data. For background electrons at the 
high energy end (~TeV) we implement the same exponential cutoff as described above. The 
results show good agreement with the observational data by e.g., ATIC and PAMELA. Fur- 
thermore, by slightly adjusting the background and a 20% decrease of the generated e + + e~ 
flux, the Fermi, H.E.S.S. and PAMELA data can also be well fitted simultaneously, as shown 
by the thin lines in Figj3j 

As a systematic check, while still assuming that pair production is the dominant in- 
teraction between CRs and photons, a constant source temperature is tried to reproduce 
the observed spectra. We find that the calculated spectra are also in agreement with the 
observational data for temperatures from 5000 K to 7000 K. 



Table 1: Parameter settings in the Monte-Carlo calculation. 
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lz 
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(10 29 cm- 2 ) 


(y r ) 




Proton 


1 x 10 4 


8.1 


0.04 


2.14 


Helium 


7 x 10 3 


12.9 


0.19 


2.02 


Iron 


2 x 10 3 


0.9 


0.58 


2.03 
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In summary we propose to explain the features of the CR spectra by the e + e~ pair 
production at the acceleration sources in an environment similar to the young SNRs. We 
show that the spectra of CRs and electrons/positrons agree well with the observations. This 
in turn provides a strong support that at least a portion of CRs might be accelerated at 
young SNRs. Our work provides a new interpretation to understand the abnormal positron 
fraction observed by PAMELA and the electrons excesses by ATIC/Fermi in the standard 
astrophysical framework. In addition, we can note in our model that only one simple set of 
the physical parameters of the CR sources seem to be sufficient to give good descriptions 
to the major properties of the data. It may indicate: 1) all the sources are "standard" and 
have similar parameters, 2) the average effect is equivalent to using single set of parameters, 
3) one single nearby source dominates the observed fluxes of CRs, and so on. If the CRs 
i ndeed come from a single ne arby source^], one would expect a sharp knee of the CR spectra 
( lErlykin fc Wolfendaldl 19971 ). It is interesting to note that t he Tibet AS7 experiment indeed 
observed a sharp knee structure (jAmenomori et al.l l2008ai ) . Furthermore, the young SNR 
acceleration model predicts a sharp cutoff of electron spectra beyond several TeV which can 
be tested by future high quality exp eriments, such as the up graded Tibet AS7 experiment 
with underground MUON detectors JAmenomori et al.ll2008bh . GRAPES^ and KASCADrfl 
Finally, we note that the current measurements of CR spectra around the knee energy are 
of great uncertainties for individual CR nuclei. If precise measurements become available in 
the future, it will be very helpful in understanding the model and in determining the model 
parameters. 



We thank Li Tipei, Ding Linkai, Lii Caidian, Wang Jianxiong, Li Haibo, Liu Siming, 
Deng Jinsong, Lu Fangjun, Zhu Shouhua, Lou Yuqing, Olivier Martineau-Huynh, Artin 
Teymourian and Wang Hanguo for helpful discussions. Hu Hong-Bo thanks Wang Wei for 
the long term help of the research work. This work is supported by the Minister of Science and 
Technology of China, Natural Sciences Foundation of China (Nos. 10725524 and 10773011), 
and by the Chinese Academy of Sciences (Nos. KJCX2-YW-N13 and K JCX3-SYW-N2) . 



2 Note, in such a case the background calculation of electrons/positrons using GALPROP might not be 
proper any more. However, phenomcnologically we can adjust the source parameters to get similar spectra 
of electrons/positrons for a single source with the continous source distribution, although there might be 
some constraints from, e.g., the diffuse gamma-rays. To build a fully self-consistent model of a single source 
scenario is beyond the scope of the present study. 

3 http: / /alpha. sci.osaka-cu.ac.jp/grapes3/index. html 

4 http://www-ik.fzk.de/KASCADE_home.html 
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Fig. 1. — Solid lines in each figure are the calculated energy spectra for protons, Helium, 
Iron and total CRs respectively, after including the pair production interactions with photons. 
The dashed lines in the upper two panels are Helium spectrum with pair production, pho- 
todisintegration and pion production, and the secondary proton spec trum due to photodisin - 



tegration of Helium. Observatio nal data a re, proton: Tibet hy brid (lAmenomori et al.l 120 06) 



Tibe t-BD dAmenomori et allboOoL JA CEE d Asakimori et al.lll998h . RUNJOB dApanasenko et al 
200lh , KASCADE dAntoni et a lihood hooj) EAS-TOP dNavarral booah. ATIC dAhnlbooah . Zat- 



JACEE (Asakimori et al 



sepin dZatsepin et al.lll993 ), SOKOL (llvanenko et al.lll993l). AMS (jAlcaraz et all 120001); Helium: 



1998), KASCADE 



2000 ). ATIC2 dPanov et al.1 12006 ) RUNJOB dApanasenko et al 



19931 ) , Ichi mura (llchimura et al 



CAPRICE (Boezio et al 



Antoni et al.l 12005), Ti bet-BD (lAmenomori et al 



200l|), SO KOL (llvanenko et al 



19931 ). BESS (lHaino et al.1 120041 ). AMS (lAguilar et al 



2003|); Iron: EAS-TOP dNavarral I2OO3J ), KASCADE dAntoni et al 



Juliusson ( Juliusson 1974 ). I chimura (Ichimura et al. 19931 ). C RN (Mueller et al 



dHarevamalll999r). RUNJOB dApanasenko et al 



2001 



2002), 



2005|), 



1991 

, SOK OL (llvanenko et al. I ll993l ): 



), Hareyama 
total: Tibet- 
Ill (|Amenomori et al.ll2008al ). Akeno ([Nagano et al.lll984l ). The normalized data are d erived by 
comb ining all data with a rescale based on the extrapolation of the direct measurements (jHorandel 



20031 ). 
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Fig. 2. — Upper-thin lines: Electron/positron spectra generated through photon-nuclei inter- 
actions for various CR components. Lower-thick lines: Propagated electron/positron spectra 
with arbitrary normalization. Also shown are the results with an exponential cu toff of the 
injection spectra with E c = 1.8, 2.1, 2.4 TeV respectively (lAharonian et al.ll2008l ). 
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Fig. 3. — Left: Total spectra for e + +e . The labels "ATIC" and "Fermi" in the 
brackets mean th e results to fit AT IC or Fer mi data respectively . Observational data 



are: CAPRICE f 
2008h . H.E. S.S. (LAharonian et al 



Boezio et al 



2000h. H EAT faarwick et all Il998h. ATIC (Ichang et al. 



Fermi-LAT flAbdo et all L2009h 



20081 ). H.E.S.S. low energy (lAharonian et al 
Right: Positron fraction. Observational data are 



2009|), 



AMS 



(Uguilar et alj 120071), HEAT9 4+95 JBarwick et al.l Il997h . HEATOO JCoutu et all boOlh . 



PAMELA f lAdriani et all 120091 ). 



